Natural variation in neuraminidase activity influences the evolutionary potential of the seasonal H1N1 lineage hemagglutinin.

The antigenic evolution of the influenza A virus hemagglutinin (HA) gene poses a major challenge for the development of vaccines capable of eliciting long-term protection. Prior efforts to understand the mechanisms that govern viral antigenic evolution mainly focus on HA in isolation, ignoring the fact that HA must act in concert with the viral neuraminidase (NA) during replication and spread. Numerous studies have demonstrated that the degree to which the receptor binding avidity of HA and receptor cleaving activity of NA are balanced with each other influences overall viral fitness. We recently showed that changes in NA activity can significantly alter the mutational fitness landscape of HA in the context of a lab-adapted virus strain. Here, we test whether natural variation in relative NA activity can influence the evolutionary potential of HA in the context of the seasonal H1N1 lineage (pdmH1N1) that has circulated in humans since the 2009 pandemic. We observed substantial variation in the relative activities of NA proteins encoded by a panel of H1N1 vaccine strains isolated between 2009 and 2019. We comprehensively assessed the effect of NA background on the HA mutational fitness landscape in the circulating pdmH1N1 lineage using deep mutational scanning and observed pronounced epistasis between NA and residues in or near the receptor binding site of HA. To determine whether NA variation could influence the antigenic evolution of HA, we performed neutralizing antibody selection experiments using a panel of monoclonal antibodies targeting different HA epitopes. We found that the specific antibody escape profiles of HA were highly contingent upon NA background. Overall, our results indicate that natural variation in NA activity plays a significant role in governing the evolutionary potential of HA in the currently circulating pdmH1N1 lineage.

Evidence of antigenic drift in the fusion machinery core of SARS-CoV-2 spike.

Antigenic drift of SARS-CoV-2 is typically defined by mutations in the N-terminal domain and receptor binding domain of spike protein. In contrast, whether antigenic drift occurs in the S2 domain remains largely elusive. Here, we perform a deep mutational scanning experiment to identify S2 mutations that affect binding of SARS-CoV-2 spike to three S2 apex public antibodies. Our results indicate that spatially diverse mutations, including D950N and Q954H, which are observed in Delta and Omicron variants, respectively, weaken the binding of spike to these antibodies. Although S2 apex antibodies are known to be nonneutralizing, we show that they confer protection in vivo through Fc-mediated effector functions. Overall, this study indicates that the S2 domain of SARS-CoV-2 spike can undergo antigenic drift, which represents a potential challenge for the development of more universal coronavirus vaccines.

Functional and antigenic characterization of SARS-CoV-2 spike fusion peptide by deep mutational scanning.

The fusion peptide of SARS-CoV-2 spike protein is functionally important for membrane fusion during virus entry and is part of a broadly neutralizing epitope. However, sequence determinants at the fusion peptide and its adjacent regions for pathogenicity and antigenicity remain elusive. In this study, we performed a series of deep mutational scanning (DMS) experiments on an S2 region spanning the fusion peptide of authentic SARS-CoV-2 in different cell lines and in the presence of broadly neutralizing antibodies. We identified mutations at residue 813 of the spike protein that reduced TMPRSS2-mediated entry with decreased virulence. In addition, we showed that an F823Y mutation, present in bat betacoronavirus HKU9 spike protein, confers resistance to broadly neutralizing antibodies. Our findings provide mechanistic insights into SARS-CoV-2 pathogenicity and also highlight a potential challenge in developing broadly protective S2-based coronavirus vaccines.

Stringent and complex sequence constraints of an IGHV1-69 broadly neutralizing antibody to influenza HA stem.

IGHV1-69 is frequently utilized by broadly neutralizing influenza antibodies to the hemagglutinin (HA) stem. These IGHV1-69 HA stem antibodies have diverse complementarity-determining region (CDR) H3 sequences. Besides, their light chains have minimal to no contact with the epitope. Consequently, sequence determinants that confer IGHV1-69 antibodies with HA stem specificity remain largely elusive. Using high-throughput experiments, this study reveals the importance of light-chain sequence for the IGHV1-69 HA stem antibody CR9114, which is the broadest influenza antibody known to date. Moreover, we demonstrate that the CDR H3 sequences from many other IGHV1-69 antibodies, including those to the HA stem, are incompatible with CR9114. Along with mutagenesis and structural analysis, our results indicate that light-chain and CDR H3 sequences coordinately determine the HA stem specificity of IGHV1-69 antibodies. Overall, this work provides molecular insights into broadly neutralizing antibody responses to influenza virus, which have important implications for universal influenza vaccine development.

Leveraging vaccination-induced protective antibodies to define conserved epitopes on influenza N2 neuraminidase.

There is growing appreciation for neuraminidase (NA) as an influenza vaccine target; however, its antigenicity remains poorly characterized. In this study, we isolated three broadly reactive N2 antibodies from the plasmablasts of a single vaccinee, including one that cross-reacts with NAs from seasonal H3N2 strains spanning five decades. Although these three antibodies have diverse germline usages, they recognize similar epitopes that are distant from the NA active site and instead involve the highly conserved underside of NA head domain. We also showed that all three antibodies confer prophylactic and therapeutic protection in vivo, due to both Fc effector functions and NA inhibition through steric hindrance. Additionally, the contribution of Fc effector functions to protection in vivo inversely correlates with viral growth inhibition activity in vitro. Overall, our findings advance the understanding of NA antibody response and provide important insights into the development of a broadly protective influenza vaccine.

Widespread impact of immunoglobulin V-gene allelic polymorphisms on antibody reactivity.

The ability of the human immune system to generate antibodies to any given antigen can be strongly influenced by immunoglobulin V-gene allelic polymorphisms. However, previous studies have provided only limited examples. Therefore, the prevalence of this phenomenon has been unclear. By analyzing >1,000 publicly available antibody-antigen structures, we show that many V-gene allelic polymorphisms in antibody paratopes are determinants for antibody binding activity. Biolayer interferometry experiments further demonstrate that paratope allelic polymorphisms on both heavy and light chains often abolish antibody binding. We also illustrate the importance of minor V-gene allelic polymorphisms with low frequency in several broadly neutralizing antibodies to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza virus. Overall, this study not only highlights the pervasive impact of V-gene allelic polymorphisms on antibody binding but also provides mechanistic insights into the variability of antibody repertoires across individuals, which in turn have important implications for vaccine development and antibody discovery.

An explainable language model for antibody specificity prediction using curated influenza hemagglutinin antibodies.

Despite decades of antibody research, it remains challenging to predict the specificity of an antibody solely based on its sequence. Two major obstacles are the lack of appropriate models and inaccessibility of datasets for model training. In this study, we curated a dataset of >5,000 influenza hemagglutinin (HA) antibodies by mining research publications and patents, which revealed many distinct sequence features between antibodies to HA head and stem domains. We then leveraged this dataset to develop a lightweight memory B cell language model (mBLM) for sequence-based antibody specificity prediction. Model explainability analysis showed that mBLM captured key sequence motifs of HA stem antibodies. Additionally, by applying mBLM to HA antibodies with unknown epitopes, we discovered and experimentally validated many HA stem antibodies. Overall, this study not only advances our molecular understanding of antibody response to influenza virus, but also provides an invaluable resource for applying deep learning to antibody research.

Stringent and complex sequence constraints of an IGHV1-69 broadly neutralizing antibody to influenza HA stem.

IGHV1-69 is frequently utilized by broadly neutralizing influenza antibodies to the hemagglutinin (HA) stem. These IGHV1-69 HA stem antibodies have diverse complementarity-determining region (CDR) H3 sequences. Besides, their light chains have minimal to no contact with the epitope. Consequently, sequence determinants that confer IGHV1-69 antibodies with HA stem specificity remain largely elusive. Using high-throughput experiments, this study revealed the importance of light chain sequence for the IGHV1-69 HA stem antibody CR9114, which is the broadest influenza antibody known to date. Moreover, we demonstrated that the CDR H3 sequences from many other IGHV1-69 antibodies, including those to HA stem, were incompatible with CR9114. Along with mutagenesis and structural analysis, our results indicate that light chain and CDR H3 sequences coordinately determine the HA stem specificity of IGHV1-69 antibodies. Overall, this work provides molecular insights into broadly neutralizing antibody responses to influenza virus, which have important implications for universal influenza vaccine development.

Widespread impact of immunoglobulin V gene allelic polymorphisms on antibody reactivity.

The ability of human immune system to generate antibodies to any given antigen can be strongly influenced by immunoglobulin V gene (IGV) allelic polymorphisms. However, previous studies have provided only a limited number of examples. Therefore, the prevalence of this phenomenon has been unclear. By analyzing >1,000 publicly available antibody-antigen structures, we show that many IGV allelic polymorphisms in antibody paratopes are determinants for antibody binding activity. Biolayer interferometry experiment further demonstrates that paratope allelic mutations on both heavy and light chain often abolish antibody binding. We also illustrate the importance of minor IGV allelic variants with low frequency in several broadly neutralizing antibodies to SARS-CoV-2 and influenza virus. Overall, this study not only highlights the pervasive impact of IGV allelic polymorphisms on antibody binding, but also provides mechanistic insights into the variability of antibody repertoires across individuals, which in turn have important implications for vaccine development and antibody discovery.

Caprin-1 binding to the critical stress granule protein G3BP1 is influenced by pH.

G3BP is the central node within stress-induced protein-RNA interaction networks known as stress granules (SGs). The SG-associated proteins Caprin-1 and USP10 bind mutually exclusively to the NTF2 domain of G3BP1, promoting and inhibiting SG formation, respectively. Herein, we present the crystal structure of G3BP1-NTF2 in complex with a Caprin-1-derived short linear motif (SLiM). Caprin-1 interacts with His-31 and His-62 within a third NTF2-binding site outside those covered by USP10, as confirmed using biochemical and biophysical-binding assays. Nano-differential scanning fluorimetry revealed reduced thermal stability of G3BP1-NTF2 at acidic pH. This destabilization was counterbalanced significantly better by bound USP10 than Caprin-1. The G3BP1/USP10 complex immunoprecipated from human U2OS cells was more resistant to acidic buffer washes than G3BP1/Caprin-1. Acidification of cellular condensates by approximately 0.5 units relative to the cytosol was detected by ratiometric fluorescence analysis of pHluorin2 fused to G3BP1. Cells expressing a Caprin-1/FGDF chimera with higher G3BP1-binding affinity had reduced Caprin-1 levels and slightly reduced condensate sizes. This unexpected finding may suggest that binding of the USP10-derived SLiM to NTF2 reduces the propensity of G3BP1 to enter condensates.

High-throughput identification of prefusion-stabilizing mutations in SARS-CoV-2 spike.

Designing prefusion-stabilized SARS-CoV-2 spike is critical for the effectiveness of COVID-19 vaccines. All COVID-19 vaccines in the US encode spike with K986P/V987P mutations to stabilize its prefusion conformation. However, contemporary methods on engineering prefusion-stabilized spike immunogens involve tedious experimental work and heavily rely on structural information. Here, we establish a systematic and unbiased method of identifying mutations that concomitantly improve expression and stabilize the prefusion conformation of the SARS-CoV-2 spike. Our method integrates a fluorescence-based fusion assay, mammalian cell display technology, and deep mutational scanning. As a proof-of-concept, we apply this method to a region in the S2 domain that includes the first heptad repeat and central helix. Our results reveal that besides K986P and V987P, several mutations simultaneously improve expression and significantly lower the fusogenicity of the spike. As prefusion stabilization is a common challenge for viral immunogen design, this work will help accelerate vaccine development against different viruses.

Mutational fitness landscape of human influenza H3N2 neuraminidase.

Influenza neuraminidase (NA) has received increasing attention as an effective vaccine target. However, its mutational tolerance is not well characterized. Here, the fitness effects of >6,000 mutations in human H3N2 NA are probed using deep mutational scanning. Our result shows that while its antigenic regions have high mutational tolerance, there are solvent-exposed regions with low mutational tolerance. We also find that protein stability is a major determinant of NA mutational fitness. The deep mutational scanning result correlates well with mutational fitness inferred from natural sequences using a protein language model, substantiating the relevance of our findings to the natural evolution of circulating strains. Additional analysis further suggests that human H3N2 NA is far from running out of mutations despite already evolving for >50 years. Overall, this study advances our understanding of the evolutionary potential of NA and the underlying biophysical constraints, which in turn provide insights into NA-based vaccine design.

Probing the biophysical constraints of SARS-CoV-2 spike N-terminal domain using deep mutational scanning.

Increasing the expression level of the SARS-CoV-2 spike (S) protein has been critical for COVID-19 vaccine development. While previous efforts largely focused on engineering the receptor-binding domain (RBD) and the S2 subunit, the amino-terminal domain (NTD) has been long overlooked because of the limited understanding of its biophysical constraints. In this study, the effects of thousands of NTD single mutations on S protein expression were quantified by deep mutational scanning. Our results revealed that in terms of S protein expression, the mutational tolerability of NTD residues was inversely correlated with their proximity to the RBD and S2. We also identified NTD mutations at the interdomain interface that increased S protein expression without altering its antigenicity. Overall, this study not only advances the understanding of the biophysical constraints of the NTD but also provides invaluable insights into S-based immunogen design.

Molecular analysis of a public cross-neutralizing antibody response to SARS-CoV-2.

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concerns (VOCs) continue to emerge, cross-neutralizing antibody responses become key toward next-generation design of a more universal COVID-19 vaccine. By analyzing published data from the literature, we report here that the combination of germline genes IGHV2-5/IGLV2-14 represents a public antibody response to the receptor-binding domain (RBD) that potently cross-neutralizes a broad range of VOCs, including Omicron and its sub-lineages. Detailed molecular analysis shows that the complementarity-determining region H3 sequences of IGHV2-5/IGLV2-14-encoded RBD antibodies have a preferred length of 11 amino acids and a conserved HxIxxI motif. In addition, these antibodies have a strong allelic preference due to an allelic polymorphism at amino acid residue 54 of IGHV2-5, which is located at the paratope. These findings have important implications for understanding cross-neutralizing antibody responses to SARS-CoV-2 and its heterogenicity at the population level as well as the development of a universal COVID-19 vaccine.

High-throughput identification of prefusion-stabilizing mutations in SARS-CoV-2 spike.

Designing prefusion-stabilized SARS-CoV-2 spike is critical for the effectiveness of COVID-19 vaccines. All COVID-19 vaccines in the US encode spike with K986P/V987P mutations to stabilize its prefusion conformation. However, contemporary methods on engineering prefusion-stabilized spike immunogens involve tedious experimental work and heavily rely on structural information. Here, we established a systematic and unbiased method of identifying mutations that concomitantly improve expression and stabilize the prefusion conformation of the SARS-CoV-2 spike. Our method integrated a fluorescence-based fusion assay, mammalian cell display technology, and deep mutational scanning. As a proof-of-concept, this method was applied to a region in the S2 domain that includes the first heptad repeat and central helix. Our results revealed that besides K986P and V987P, several mutations simultaneously improved expression and significantly lowered the fusogenicity of the spike. As prefusion stabilization is a common challenge for viral immunogen design, this work will help accelerate vaccine development against different viruses.

Prevalence and mechanisms of evolutionary contingency in human influenza H3N2 neuraminidase.

Neuraminidase (NA) of human influenza H3N2 virus has evolved rapidly and been accumulating mutations for more than half-century. However, biophysical constraints that govern the evolutionary trajectories of NA remain largely elusive. Here, we show that among 70 natural mutations that are present in the NA of a recent human H3N2 strain, >10% are deleterious for an ancestral strain. By mapping the permissive mutations using combinatorial mutagenesis and next-generation sequencing, an extensive epistatic network is revealed. Biophysical and structural analyses further demonstrate that certain epistatic interactions can be explained by non-additive stability effect, which in turn modulates membrane trafficking and enzymatic activity of NA. Additionally, our results suggest that other biophysical mechanisms also contribute to epistasis in NA evolution. Overall, these findings not only provide mechanistic insights into the evolution of human influenza NA and elucidate its sequence-structure-function relationship, but also have important implications for the development of next-generation influenza vaccines.

Egg-adaptive mutations of human influenza H3N2 virus are contingent on natural evolution.

Egg-adaptive mutations in influenza hemagglutinin (HA) often emerge during the production of egg-based seasonal influenza vaccines, which contribute to the largest share in the global influenza vaccine market. While some egg-adaptive mutations have minimal impact on the HA antigenicity (e.g. G186V), others can alter it (e.g. L194P). Here, we show that the preference of egg-adaptive mutation in human H3N2 HA is strain-dependent. In particular, Thr160 and Asn190, which are found in many recent H3N2 strains, restrict the emergence of L194P but not G186V. Our results further suggest that natural amino acid variants at other HA residues also play a role in determining the preference of egg-adaptive mutation. Consistently, recent human H3N2 strains from different clades acquire different mutations during egg passaging. Overall, these results demonstrate that natural mutations in human H3N2 HA can influence the preference of egg-adaptation mutation, which has important implications in seed strain selection for egg-based influenza vaccine.

The evolutionary potential of influenza A virus hemagglutinin is highly constrained by epistatic interactions with neuraminidase.

Antigenic evolution of the influenza A virus (IAV) hemagglutinin (HA) gene limits efforts to effectively control the spread of the virus in the population. Efforts to understand the mechanisms governing HA antigenic evolution typically examine the HA gene in isolation. This can ignore the importance of balancing HA receptor binding activities with the receptor-destroying activities of the viral neuraminidase (NA) to maintain viral fitness. We hypothesize that the need to maintain functional balance with NA significantly constrains the evolutionary potential of the HA. We use deep mutational scanning and show that variation in NA activity significantly reshapes the HA fitness landscape by modulating the overall mutational robustness of HA. Consistent with this, we observe that different NA backgrounds support the emergence of distinct repertoires of HA escape variants under neutralizing antibody pressure. Our results reveal a critical role for intersegment epistasis in influencing the evolutionary potential of the HA gene.

Interactions between Influenza A Virus Nucleoprotein and Gene Segment Untranslated Regions Facilitate Selective Modulation of Viral Gene Expression.

The influenza A virus (IAV) genome is divided into eight negative-sense, single-stranded RNA segments. Each segment exhibits a unique level and temporal pattern of expression; however, the exact mechanisms underlying the patterns of individual gene segment expression are poorly understood. We previously demonstrated that a single substitution in the viral nucleoprotein (NP:F346S) selectively modulates neuraminidase (NA) gene segment expression while leaving other segments largely unaffected. Given what is currently known about NP function, there is no obvious explanation for how changes in NP can selectively modulate the replication of individual gene segments. In this study, we found that the specificity of this effect for the NA segment is virus strain specific and depends on the untranslated region (UTR) sequences of the NA segment. While the NP:F346S substitution did not significantly alter the RNA binding or oligomerization activities of NP in vitro, it specifically decreased the ability of NP to promote NA segment viral RNA (vRNA) synthesis. In addition to NP residue F346, we identified two other adjacent aromatic residues in NP (Y385 and F479) capable of similarly regulating NA gene segment expression, suggesting a larger role for this domain in gene-segment specific regulation. Our findings reveal a novel role for NP in selective regulation of viral gene segment replication and provide a framework for understanding how the expression patterns of individual viral gene segments can be modulated during adaptation to new host environments. IMPORTANCE Influenza A virus (IAV) is a respiratory pathogen that remains a significant source of morbidity and mortality. Escape from host immunity or emergence into new host species often requires mutations that modulate the functional activities of the viral glycoproteins hemagglutinin (HA) and neuraminidase (NA), which are responsible for virus attachment to and release from host cells, respectively. Maintaining the functional balance between the activities of HA and NA is required for fitness across multiple host systems. Thus, selective modulation of viral gene expression patterns may be a key determinant of viral immune escape and cross-species transmission potential. We identified a novel mechanism by which the viral nucleoprotein (NP) gene can selectively modulate NA segment replication and gene expression through interactions with the segment UTRs. Our work highlights an unexpected role for NP in selective regulation of expression from the individual IAV gene segments.